Flow of He II due to an Oscillating Grid in the Low Temperature Limit 
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The macroscopic flow properties of pure He II are probed in the limit of zero temperature using an 
oscillating grid. With increasing oscillation amplitude the initially pure superflow changes abruptly; 
the resonant frequency decreases and the response becomes strongly nonlinear, attributable to the 
growth of a boundary layer of quantized vortices that increases the effective mass of the grid. On 
further increase of oscillation amplitude, the flow undergoes a transition to turbulence. 

PACS numbers: , 47.37.+q, 47.27.Cn, 47.15. Cb, 67.40.Vs 



Although the exotic flow properties of He II have been 
a subject of intensive investigation ever since the discov- 
ery of superfluidity, accumulating a vast amount of ex- 
perimental data and theoretical knowledge many 
important features still remain to be explained. In the 
limit of low velocity. He II flow is very well described in 
terms of Landau's two-fluid model. On increasing the 
flow velocity beyond a certain threshold, however, quan- 
tized vortices appear. Their presence couples the orig- 
inally independent normal and superfluid velocity fields 
via a mutual friction. Above ^1.2 K, where He II con- 
tains an appreciable proportion of normal fluid, numer- 
ous investigators observed that on exceeding a suitably 
defined Reynolds number. He II flow acquires an increas- 
ingly classical character: (i) the surface of a rotating 
bucket forms a nearly classical meniscus Q ; (ii) flow past 
a microsphere displays both laminar and turbulent drag 
, as well as (iii) a drag crisis ^ ; (iv) the energy spec- 
trum of turbulent He II involves an inertial range with 
a classical K41 Kolmogorov roll-off exponent of —5/3 
; (v) the decay of both grid-generated and counterflow 
turbulence displays a classical character 0, 0, 0, 
Although these features are seen over a wide range of 
normal fluid/superfluid density ratios, it is impossible to 
exclude the possibility that this classical-like behaviour 
is caused by the presence of the viscous normal fluid. 
Following oscillating experiments with a tiny sphere Q 
and a very thin vibrating wire j^, yielding interesting 
results attributable to single vortices, there is a clear call 
to study the macroscopic flow properties of He II in the 
T ^ limit where normal fluid is (almost) absent and 
pure superflow can be investigated systematically. 

Our tool for performing such investigations is the thin 
electroformed nickel grid |T3| of 70% transparency and 
mesh size of 127 /im shown schematically in Fig. ^ It 
constitutes a circular membrane 2R = 8 cm in diameter, 
tightly stretched on a circular mild steel carrier, mounted 
horizontally, equidistant from two gold-plated plane cop- 
per electrodes, each drilled with 170 holes of 2 mm ID 
to connect the grid space to the rest of the experimental 
chamber. The electrodes are separated by = 2 mm. 




FIG. 1: Schematic drawing showing the geometry of the grid 
and an electron microscope picture of the actual grid (sample 
cut off the same sheet). 



immersed in a 1.5 liter volume of isotopically pure ^He 
in an experimental cell attached to the mixing chamber 
of the dilution refrigerator. The apparatus can be op- 
erated at all pressures up to the solidification pressure. 
A static potential, typically Vq = 500 volts, is applied 
to the grid. A driving potential Vi = Vwcosujt [Vw 
^ Vb) applied to the upper electrode produces a driving 
force on the grid of form fd — 2eoer'n' R^VVi/cP, where 
Eq and Er denote respectively the relative permittivities 
of free space and the liquid '''He. The grid can thus be 
considered as an oscillating membrane |l3| under uni- 
form tension. In this study we approximate its motion 
as one-dimensional and assume that the oscillation ampli- 
tude is uniform across its area. Oscillations of amplitude 
AD induce a signal of amplitude V2 = VAD/{2d) on 
the lower electrode, which can be monitored with a lock- 
in amplifier (which we use to investigate the low-drive 
linear response of the grid) or directly with a memory 
oscilloscope, allowing visualization of transient processes 
resulting from the strongly nonlinear response at higher 
drives. Allowing for a reduction in the induced voltage V2 
by a factor of (1 + Cc/C)-^ « 0.065, where = 700 pF 
is the capacitance of the connecting cable and C = 47 pF 
is the capacitance between the grid and the lower copper 
electrode, the response amplitude IV2I provides a direct 
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FIG. 2: Resonance curves measured at 10 bar and nominally 50 inK using the memory oscilloscope for various drive levels (in 
Vpp) as indicated. Doubled symbols in some places indicate where stable beatings were observed while carefully sweeping the 
drive frequency, typically in steps of 0.001 Hz. The solid line represents the resonant response to a 0.05 Vpp drive in vacuum 
(measured using the lock-in amplifier, upper frequency axes). 



measure of the amplitude and spatially averaged peak 
velocity \vg\ = \ujAD\ of the oscillating grid. 

At the lowest oscillation amplitudes, the grid displays 
linear behaviour in that its frequency response to the 
drive is a Lorcntzian of narrow width caused predom- 
inantly by nuisance damping, with a quality Q factor 
typically exceeding 5000. The resonant frequency /res 
can be altered temporarily by increasing the pressure, 
which we attribute to quantized vortices generated e.g. 
by the jet from the filling capillary: they probably be- 
come pinned to inhomogeneities on the grid surface, and 
between the grid and surrounding electrodes. Moving the 
grid violently at high drive amplitude, however, is found 
to stabilize /res typically to within ±0.1 Hz, presumably 
by shaking off most of this pinned vorticity. 

Fig. 121 shows our central experimental observation, 
which for the "cleaned" grid has been repeated for several 
pressures spanning 0.3 < p < 24.8 bar without apprecia- 
ble change. As the drive increases and the grid amplitude 
reaches the first threshold (typically 10-20 mVpp, corre- 
sponding to a mean grid velocity of 0.3 < Wg^"* < 0.6 
cm/s), the oscillation amplitude at resonance continues 



to rise in proportion to the drive (Fig. |2Jl, but the res- 
onant frequency suddenly starts decreasing and the res- 
onance curves acquire highly nonlinear features. If one 
follows e.g. the resonance curve for the 0.05 Vpp drive 
down from 1091 Hz (Fig.[21) by slowly sweeping the drive 
frequency (in practice, digitally, in steps of 0.001 Hz), 
the system displays the usual nearly Lorentzian station- 
ary response while the amplitude remains below the first 
threshold. On further decreasing the drive frequency, the 
response amplitude increases above this first threshold 
until, shortly below 1086 Hz the amplitude suddenly col- 
lapses down to a lower stable branch. On increasing the 
drive frequency again from this point the system stays on 
the lower branch until about 1087.15 Hz, where a transi- 
tion to the stable upper branch occurs. These hysteretic 
loops are robust to temperature increase, at least up to 
our maximum of 130 niK. 

Within the hysteretic parameter range, beat phenom- 
ena are occasionally seen between two apparently sta- 
ble amplitude levels, with envelopes of period ~1 s (see, 
eg.. Fig. 13 drive 0.05 Vpp around 1086.5 Hz); once es- 
tablished, they are stable on the scale of hours. Small 
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FIG. 3: Response amplitude of the grid versus the drive level 
at resonance (main figure) and at fixed frequencies (inset). 
For explanation of hysteresis, see text. The full lines indicate 
the linear and square-root responses. Note that the main 
figure and inset refer to different pressures. 



changes of frequency do not kill them, but modify the 
upper and lower amplitude levels between which beating 
occurs. Only with further change of the driving frequency 
does the beating disappear. On restoring the original 
frequency the response is not completely reproducible; 
within some small frequency range beating might not ap- 
pear at all; and sometimes it reappeared only later, or not 
at all. Despite considerable effort, we did not succeed in 
establishing any fully repeatable pattern or well-defined 
conditions for the appearance of beats . 

The fall in frequency with increasing drive (Fig. ^ 
reaches typically 2 Hz for all investigated pressures (see 
Fig. ceasing at a second threshold amplitude (typ- 
ically 200 mVpp corresponding to a mean grid velocity 

of about Vg ' X 6 cm/s). Observation of the two well- 
defined resonant frequencies, at all investigated pressures 
shifted by « 2 Hz, is a remarkable feature of the superflow 
that, to our knowledge, has not previously been reported. 
With further increase in drive, the oscillation amplitude 
at resonance initially remains almost constant, while the 
widths of the resonance curves increase rapidly (see Fig. 
121). Only for drive levels exceeding by about an order of 
magnitude that for the second threshold does the ampli- 
tude at resonance grow again; this time approximately in 
proportion to the square-root of the drive, as shown in 
Fig. El In this high-drive regime the linewidth increases 
rapidly while the resonant frequency decreases gradu- 
ally, qualitatively in the manner expected for increasing 
damping. Pronounced beatings were never observed in 
this high-drive regime, ff the oscillation amplitude is 
measured as a function of drive level at fixed frequency, 
one observes clear hysteresis within the frequency range 
containing the two stable branches of the grid response 
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FIG. 4: Frequency as a function of He II density at low drive 
levels (circles) and for the second critical threshold (trian- 
gles) The straight line extrapolates through the zero-density 
(vacuum) resonant frequency of /o=1117.2 Hz. 



(Fig. O . Outside this range, the drive dependences are 
nonlinear, but without any hysteresis. 

Measurements of the grid response at low temperature 
in vacuum (smooth curve in Fig.j^J demonstrate that the 
presence of the He II has two effects. First, it down-shifts 
the low-drive resonant frequency by about 30 Hz from 
its vacuum value of /o = (1117.2 ± 0.05) Hz, depending 
on the pressure p in the cell (see Fig. 0}. This classi- 
cal effect arises from hydrodynamic enhancement of the 
mass M of the grid by AMh = l3gHe{p)M / qg, [HIi^. 
The measured /o and low drive resonant frequencies yield 
(3 = 3.01 ± 0.05. Secondly, it is indeed the He II that is 
responsible for the observed nonlinearities |l(il |. 

We believe that the behaviour of the oscillating grid 
ought to be understandable in terms of quantized vor- 
tices generated in He II by its motion. As one possibility, 
we may speculate that, on exceeding the first threshold, 
a "bound ary layer" of vortex loops (perhaps of horseshoe 
form [13,113) grows on the grid, thus increasing its effec- 
tive mass. Correspondingly, the resonant frequency shifts 
down and the resonance curves acquire strongly nonlinear 
features. Provided the vortex loops remain pinned and 
that they do not reconnect to create free vorticity, the 
mass enhancement would be non-dissipative, as observed 
(see below). We note that small mass enhancements were 
observed earlier for a sphere vibrating in He II above 1 K 
[l^. but without associated critical thresholds. 

It is interesting to characterize the first threshold by 
a superfluid Reynolds number Rcs = C/maxG/K, where 
C^max stands for a peak critical flow velocity through the 
grid window j23|, G is its linear size and k denotes the 
circulation quantum. Observed values of Rcg 10 com- 
pare well with the critical Rcs = UD/k ~ 20 {U \s the 
transport velocity and D the diameter of the pipe) found 
as a temperature independent threshold in pipe flow of 
He II at much higher temperature when the flow of the 
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normal component was inhibited by superleaks at both 
ends of the pipe . The generation of quantized vortex 
loops on the surface of the grid by macroscopic super- 
flow around it p robably involves growth from remanent 
vortex lines |22| , or from a "plasma" of half vortex rings 
[l7|. This idea is supported by the following observa- 
tion. While measuring the response dependence on in- 
creasing the drive amplitude at the resonant frequency, 
the system encounters a nucleation problem when pass- 
ing the first threshold: on some occasions the resonance 
response stopped growing with increasing the drive level, 
and jumped on the response/drive curve (Fig. |3J only 
later. With decreasing drives this feature disappears, and 
the response remains proportional to the drive level. 

The frequency down-shift Af{p) ~ /ics — fsh ^ 2 Hz 
between the two critical thresholds might be considered 
in terms of a boundary layer of thickness A that en- 
hances the hydrodynamic effective mass of the grid by 
AMg = AXgHeip), where A denotes the surface area 
of the grid. Requiring that the downward shift of the 
resonance frequency corresponds to those observed ex- 
perimentally leads to 



A = -. ^ 



0.53 ±0.05 ^im 



which would correspond to an effective kinematic viscos- 
ity i^eff ~ 10~^ cm^/s. Extending the analogy to classical 
fluids, we may estimate the expected resonant linewidth 
due to the drag of such a hypothetical viscous fluid. A 
straightforward calculation |23| based on estimating the 
flow velocity gradient in the direction inside the fluid by 
{vg)/X leads to a linewidth ^1 Hz, consistent with res- 
onances for (e.g. Fig. [3 0.2 and 0.5 Vpp data) near the 
second critical threshold. However, this picture seems 
unable to account for the re-entrant character of the non- 
linear resonances because the resultant viscous broaden- 
ing would tend to smear out the effect in question. 

We infer that, beyond the second critical threshold 
(Fig.[2Jl, the vortex loops start to reconnect. Free vortic- 
ity is then shed into the liquid, corresponding to dissipa- 
tion of the vibrational energy and leading to the observed 
increase in linewidth. Under these flow conditions. He II 
behaves in close analogy with a classical Navier-Stokes 
fluid: the square-root behaviour (Fig. ^ of the resonant 
response as a function of drive amplitude in this range 
is typical of classical turbulent drag scaling, and it is 
therefore likely that this threshold marks the onset of 
turbulence. After leaving t he g rid, the quantized vor- 
ticity probably "evaporates" [2J| as recently proposed in 
connection with turbulence in superfluid '^He-B |25| . 

The particular challenge posed by these results is to de- 
velop a quantitative description of vortex dynamics in the 
macroscopic flow around the moving grid, showing how it 



can give rise to the observed amplitude-dependent mass 
enhancements and re-entrant resonance curves, without 
dissipation. But the problem also carries wider signifi- 
cance, and will repay careful study, because the dynam- 
ical processes in question are probably fundamental to 
the generation of quantum turbulence. 
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